Decoherence and dephasing in multilevel systems interacting with thermal 

environment 
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We examine the effect of multilevels on decoherence and dephasing properties of a quantum system 
consisting of a non-ideal two level subspace, identified as the qubit and a finite set of higher energy 
levels above this qubit subspace. The whole system is under interaction with an environmental bath 
through a Caldeira-Leggett type coupling. The model interaction we use can generate nonnegligible 
couplings between the qubit states and the higher levels upto N ^ 10. In contrast to the pure two- 
level system, in a multilevel system the quantum information may leak out of the qubit subspace 
through nonresonant as well as resonant excitations induced by the environment. The decoherence 
properties of the qubit subspace is examined numerically using the master equation formalism of the 
system's reduced density matrix. We numerically examine the relaxation and dephasing times as 
the environmental frequency spectrum, the environmental temperature, and the multilevel system 
parameters are varied. We observe the influence of all energy scales in the noise spectrum on 
the short time dynamics implying the dominance of nonresonant transitions at short times. The 
relaxation and dephasing times calculated, strongly depend on A'^ for 4 < A'^ < 10 and saturate for 
10 < N . We also examine double degenerate systems with 4 < A'^ and observe a strong suppression 
(almost by two orders of magnitude) of the low temperature relaxation and dephasing rates. 

An important observation for 4 < A'^ in doubly degenerate energy configuration is that, we find a 
strong suppression of the RD rates for such systems relative to the singly degenerate ones. These 
results are also compared qualitatively with the relaxation rates found from the Fermi Golden Rule. 

PACS numbers: 03.65.X,85.25.C,85.25.D 



I. INTRODUCTION 

Currently a large number of model approaches are 
present for formulating the decoherence phenomena in 
the literature. The original Caldeira-Leggett model^ is 
based on a quantum system under the influence of a dou- 
ble well tunneling potential with a linear coupling to an 
infinite bath of harmonic oscillators. If the potential is 
sufficiently smooth and the high energy levels are suffi- 
ciently above the tunneling barrier this original model 
is normally represented as a two level system^ (2LS) in- 
teracting with a bosonic environmental bath (spin-boson 
model) . An incomplete list of this wide literature is pro- 
vided inSi^. Another popular model of decoherence is 
the central spin system in which central 2LS couples to 
large number of environmental two level systems. The 
pros and cons of these two rival models have been exten- 
sively studiec^.. 

Realistically, and aside from the genuine 2LS, a large 
majority of physical systems suggested for qubit is far 
from being ideal. In a quantum computational environ- 
ment, the parameters of the physical systems are manip- 
ulated to perform the gate operations. For instance, in 
a multilevel system (MLS), short time pulses used in the 
manipulation of the states in the qubit subspace induce 
nonresonant transitions out of the qubit subspace. That 
nonresonant transitions contribute to the decoherence of 
the quantum system was recently addressed by Tian and 
Lloyd^. They suggest that after these transitions are in- 
duced an optimized sequence of controlled pulses can be 
applied to cancel the nonresonant effects at arbitrary pre- 
cision. The idea being physically correct, requires an ad- 



ditional fine knitting of error correction which undoubt- 
edly costs computational time. On the other hand, one 
may address the same issue by seeking for an alternative 
solution: Can one understand the effect of the higher lev- 
els on decoherence in a well-parameterized MLS coupled 
to an environment? 

The MLS can itself be manifestly A^-levelled or a trun- 
cated approximation of a larger system with much higher 
number of levels. Well known examples of both cases 
have been known. For the former, organic molecules 
with certain discrete rotational symmetries and ground 
state low temperature configurations of single polimer- 
ized chains are good examples. The vibrational energy 
spectra of atoms and molecules is a good example for the 
latter. 

These type of realistic MLS can be found for instance 
in already well-examined superconducting systems such 
as the rf-SQUID in the charge, flux or phase regimes. 
We remark however, that a concise treatment of the de- 
coherence effects in MLS has not been fully developed. 
This work is planned to be a modest step forward in that 
direction. 

In section II we give an introduction of the model MLS 
used in the present work. Here we merely concentrate 
on the properties of the environmentally induced dipole 
matrix elements. Section III recalls the reduced density 
matrix (RDM) master equation formalism and adapts it 
for the coupling of the MLS to the environment. The 
noise correlator, which is considered to be in thermal 
equilibrium, and the system-noise kernel, for which no 
Markovian assumption is made, are defined in section 
III. A. The results are presented together with the ear- 
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lier observations of the 2LS (in section III.B) to allow a 
comparison with some of the established facts. The MLS 
with three or higher levels are examined in section III.C 
separately for = 3, iV = 4 and 4 < iV. The section 
III.C and the following section IV comprise most of the 
original results of the manuscript. In section IV, the re- 
laxation times for MLS, which we produced numerically 
in section III, are compared with the estimates found by 
using the semianalytic Fermi Golden Rule (FGR). 

II. THE MULTILEVEL SYSTEM 

In principle the majority of the well established meth- 
ods (particularly the influence functional) used in the 
literature automatically accomodate multilevel dynam- 
ics. The results are generally hoped to be true for 2LS in 
the WKB limit at sufficiently low temperatures^'^. Such 
techniques are also often preferable since they allow ex- 
plicit analytic expressions for the decoherence times as 
functions of the system's parameters. On the other hand, 
exact methods are also available on pure 2LS'*. However, 
we need to accomodate in our parameters explicit de- 
generacy factors as well as the Hilbert space dimensions. 
For the latter, the influence functional formalism has no 
parametric dependence. In this work, in order to retain in 
our calculations the dependence on the system's Hilbert 
space dimension, we resort to the system's eigenenergy 
basis representation. 

Our MLS model is an rf-SQUID operating in the quan- 
tum coherence regime given by the Hamiltonian 

Hs/ihno) ^^[-dl + {z~ipM,sf]+l3 cos{-fz) (1) 

where Hq = 271 /\/LC is the harmonic energy with L be- 
ing the inductance of the SQUID loop and C is the effec- 
tive capacitance of the Josephson junction, (3 = Ej/MIq 
is the dimensionless ratio of the Josephson energy Ej 

to the harmonic energy, 7 — ^^/^(f^) is a dimen- 
sionless scale parameter, ipbias — '^'Tr^uas/'^o is the ef- 
fective bias in the flux which is applicable in a current 
bias junction, and z — 27r$/(7<I>o) is the flux ($) de- 
pendent dimensionless dynamical variable ($0 = hc/2e 
is the superconducting flux quantum). Truncating the 
Hilbert space dimension at N in the energy eigenbasis, 
the system Hamiltonian in is written in the diagonal 
form 

N-l 

ns=Y. Er.i{C})m,n){{C},n\ (2) 

n=0 

where {(,} describes the set of system parameters 
rio, f3,j,ipbias where £'„({(}), and |{C},"-) are respec- 
tively the parameter dependent eigenenergies and eigen- 
vectors of the MLS. This set of parameters is sufficiently 
general to accomodate for all possible effects includ- 
ing the degeneracy in the qubit subspace, the symme- 
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FIG. 1: The double- well potential and the eigenenergy con- 
figurations corresponding to the singly degenerate (SD) case. 
Here the numerical values of the dimensionless parameters for 
this SD configuration are /3 ~ 1.616 and 7 ~ 1.753. The har- 
monic energy scale and the degeneracy factors are respectively 
hQ,o — 10~^eV and r]2 = 10^. Here the numerical values of 
the dimensionless parameters are /? ~ 1.616 and 7 ~ 1.753. 



try of the wavefunctions etc. [we define the degener- 
acy factor by 77 = {E2 - Ei)/{Ei - Eo) for MLS and 
= {El + Eq)/{Ei - Eo) for 2LS]. These three pa- 
rameters VLq, (3,ipbias control the high energy harmonic 
spectrum, the low energy anharmonic spectrum and the 
refiection symmetry of the rf-SQUID potential respec- 
tively. At low energies, a simple numerical diagonaliza- 
tion reveals that there are low lying eigenenergy con- 
figurations within the double well regime in which the 
SQUID potential is strongly anharmonic. The param- 
eters of the potential can therefore be manipulated to 
search within this regime for those configurations satis- 
fying optimal qubit conditions. An interesting case here 
is to find highly degenerated levels corresponding to the 
first two eigenstates for the symmetric double well po- 
tential (i.e. iphias = 0). This particular case has been 
extensively examined previously for 2LS using semiclas- 
sical methods with an arbitrarily weak tunneling between 
the wellsS'^. Another configuration that turns out to be 
important in our calculations is the doubly degenerate 
(DD) configuration for systems with 4 < iV in which 
the first four levels are pairwise degenerate with large 
degeneracy factors. The double-well potential and ener- 
gies corresponding to both SD and DD configurations are 
shown in Fig's 1^ and (0) respectively. 

We numerically find that, the relaxation/dephasing 
(RD) times for the MLS can be controlled by the degener- 
acy rj. Normally, degeneracy is also crucial in controlling 
the dynamical tunneling rates. In our calculations how- 
ever we directly use the system eigenstates. Therefore for 
the highly degenerate configurations bare tunneling be- 
tween symmetric and anti-symmetric parts of the wave- 
functions can be neglected to a large extend. This turns 
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FIG. 2: The double-well potential and the eigenenergy con- 
figurations corresponding to the doubly degenerate (DD) case 
considered throughout the manuscript. The numerical values 
of the dimensionless parameters corresponding to this DD 
configuration are /3 ~ 0.772 and 7 ~ 2.187. The harmonic 
energy scale and the degeneracy factors are is hflo — 10^^ 
eV, 772 ~ 10'', r? ~ 3 X 10^ 



out to be especially important for quantum computation 
in the sense that once the computation is finished the 
wavefunction can be maximally localized in one of the 
double wells before any measurement or read-out pro- 
cess. 

Although we use the truncated rf-SQUID as the N-level 
model quantum system to study decoherence effects, our 
treatment is not at the microscopic level. The rf-SQUID 
is shown to be an ideal model to study multilevel effects 
due to the fact that, the transitional dipole couplings be- 
tween the low lying energy states and the high levels are 
nonncgligiblc [sec Fig.©]. Any other physical Hamilto- 
nian with sufficient number of adjustable parameters as 
well as nonnegligible dipole couplings would be suitable 
for the calculations presented here. 

In the rest of the paper the harmonic energy fig = 
2tt/\^LC is a free parameter which is used for scaling 
energy and time. 



A. Coupling to Noise 

The system-noise interaction is considered to be a 
Caldeira-Leggett type inductive coupling between the 
squid's macroscopic flux coordinate z and the envi- 
ronmental flux-like coordinate ipe expanded in harmonic 
environmental modes as (pe = J2k (^Lfe + ^fe)- The sys- 
tem noise interaction is simply Hint — ^ zipe where a is 
some number representing the strength of the inductive 
coupling {a is to be normalized by Qq for a dimension- 
less coupling). In the diagonal basis |{C}, of the model 




n (truncated Hilbert space dimension) 



FIG. 3: A few non-zero dipole matrix elements (2:)n,m of the 
coupling of the rf-SQUID to a fiux noise versus the truncated 
Hilbert space dimension [calculated in the eigenenergy basis 
of iQ] 

system the interaction Hamiltonian is given by 

Wmt = f E W-l4C})(r-{C}|^e , (3) 
r,s— 

Here {z)rs = ({C} ''l'Z|{C} s) Siie the noise induced pertur- 
bative dipole matrix elements of the macroscopic system 
coordinate z in the MLS's diagonal basis in For the 
model MLS described by , the dipole matrix elements 
are real and symmetric. 

The rf-SQUID poses a general example in which the 
multilevelledness of the master system manifests itself by 
finite dipole transition matrix elements for both the sym- 
metric (i.e. ifibias — 0) and asymmetric (i.e. (ptias 7^ 0) 
potential configurations. In the symmetric contribution 
even parity transitions vanish which results in manifestly 
off-diagonal system-noise coupling. Physically, this is in 
contrast to the most popular models used in the litera- 
ture. On the other hand, when the potential is tilted, the 
parity selection no more holds by which finite diagonal 
couplings are also created. 

In Figl^lthe noise induced couplings for an asymmet- 
ric potential is plotted as a function of the truncated 
Hilbert space dimension. Data indicate that the induced 
dipole strengths between the qubit and the higher energy 
states are comparable to those among the qubit states. 
Therefore, the high energy transitions cannot be triv- 
ially ignored. The high energy transitions normally ap- 
pear as a result of resonant interactions with the high 
energy sector of the noise spectrum under long interac- 
tion times. However, in the reduced system these transi- 
tions appear in the short time dynamics as well, and the 
short time dynamics is dominated by the nonresonant 
processes. Considering that decoherence is dominantly 
affected by the short time behaviour, the nonresonant 
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processes are expected to have observable effects in the 
decoherence properties of the RDM. Indeed we observe 
these effects in the solution of the master equation for 
the MLS (section III.B and C). 

The next is to consider the environmental spectrum 
and the availability of the bath frequencies for these exci- 
tations. Regarding this, we consider a thermal Gaussian 
environment spectrum 

I{lo) = Lo^+'' e-'^'/'*^' coth(w/2T) (4) 

where A is the effective noise cutoff frequency and v 
describes the subohmic (i.e. < 0), superohmic (i.e. 
v > 0), and ohmic (i.e. v = 0) character of the spec- 
trum. The three environmental parameters A, T deter- 
mine the sectors of the spectrum where the system-noise 
coupling is most effective. For —l^v (extreme sub- 
ohmic), two regions are of particular importance: a) at 
sufficiently low temperatures and high cutoff correspond- 
ing to T ^ w < A, the dominant mechanism of relaxation 
is through spontaneous deexcitations^. We call this re- 
gion region-1] b) at high temperatures and high cutoff 
the region < w <C min(A, T) provides a wider range of 
strong environmental couplings which we call as region- 

II. If the character of the spectrum is more like ohmic 
or superohmic, i.e. ~ or ~ 1 respectively, there 
is lesser room for dexcitations as the availability of the 
low energy modes is suppressed. Therefore, in the ohmic 
and superohmic regimes, the region-11 dominates the RD 
phenomena. 

Another feature of ^ is related to the majority of 
critical crossover behaviour in the vicinity of = as 
predicted earlier by Leggett et al. and depicted in FigQ] 
In this figure, the ~ is a critical vicinity in the Ohmic 
region separating the subohmic < regime from that 
< i^. In the subohmic regime I{uj)/ {2Ty~^'^ is very 
small except for vanishingly small frequencies (a;/2T <ti 
1). Whereas, in the regime < v the maximum value of 
is observed at higher frequencies uj ~ 2A^ (1 + i^) /2 
with an intensity proportional to (A/T)^+'^. 

III. MASTER EQUATION AND THE REDUCED 

DENSITY MATRIX FOR THE MLS 

In the study of decoherence effects due to the weak 
environmental influences, one conventional way is to cal- 
culate the time dependent RDM elements by solving 
the master equation. This formalism has been known 
since the independent works of Bloch, Redfield and Fano 
(BRF)2i on spin magnetic resonance and widely applied 
to the spin-boson systems for which many standard ref- 
erences exisfeiS. The standard BRF formalism assumes 
Markov conditions for the solution of the master equa- 
tions, which leads to exactly solvable results for 2LSii. 
However, the Markovian assumption is not free of draw- 
backs and that was questioned originally ir.^^ and lately 
in^ as well as'^- in the context of spin magnetic resonance 
and relaxation. 




FIG. 4: Tlie variation of I{uj) in @ versus cj and f for 
-1 < ly < 1 parameterized by A/T = 10,50,100 (from the 
innermost to the outermost surfaces respectively). 

In this work, the system noise kernel is treated with 
its most general non-Markovian character—. The time 
evolution of the RDM is obtained in the interaction rep- 
resentation by 

-ihj^Xt)^[pit),nrntm (5) 

where ^ denotes the interaction picture. In the context of 
decoherence, we will give more emphasis on short obser- 
vational times in the solution of (O . A convenient way to 
proceed is then to apply the Born approximation in which 
the full density matrix is initially a product of the system 

and environmental ones (i.e. p(0) = '(0)(^p' '(0))and 
at any later and short time approximately separates as 

The exact iterative solution of (0) including the second 
order in the interaction with the partial trace performed 
over the environmental degrees of freedom yields for the 
RDM the integro-differential equation 

j/nlit) = - f dt' J2 KrMpi'jit') (6) 

in which we adopt the model interaction Hamiltonian Q 
for the system-noise kernel which is found to be 

— ^|-^(^ ~ t')[{ztZt')nr5s,m ~ {zt')nr{zt) s m] (J) 
+J-*{t — t')[{zt'Zt)ms5r,n — {zt)nr{zt') sm]^ 

Note the the kernel depends on two times as a signature 
of the non Markovian treatment and it is not time trans- 
lationally invariant. Here J-{t — t') = T*{t' — t) is the 
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complex noise correlation function 



and 



N-l 

k,t=0 



(8) 



(9) 



is the time dependent dipole operator in the interaction 
picture where Ek^ \ k) comprise the eigensolution of the 
model system. Expanding the noise field if>e in the inde- 
pendent harmonic modes and calculating (jSJ in thermal 
equilibrium one obtains the standard thermal noise cor- 
relator 

T{t-t') = 2^77^,[coth(a)A./2T) cofiuok{t-t')-isi\iuok{t-t')] 

k 

(10) 

The Markovian versus non Markovian character of the so- 
lution of lO is determined in the weak system-noise inter- 
action limit by the competition of three time scales: tb, 
noise correlation time scale, and Tdep, the relaxation 
and dephasing time scales^^ of the reduced system respec- 
tively. The noise correlation time scale is found roughly 
from by the thermal Gaussian bath spectral width as 
tb — 1/A. At the Markovian limit, the environmental 
correlation time tb — A^^ is much smaller than the RD 
times. For two level systems this condition can be met 
provided that the system noise-coupling is sufficiently 
small. However, for MLS, the question of whether the 
Markovian condition is satisfied is more nontrivial. The 
basic reason is that in the MLS there is a larger number 
of time scales and decay channels of which presence may 
considerably reduce the effective decoherence times. 

In this work, the numerical solution of © is performed 
by discretizing time in steps At = 10~^A^^. The Her- 
miticity and the normalization of the RDM at each time 
step is maintained within an accuracy of 10~^^. 



A. The system-noise Kernel 

In the model Hamiltonian defined in ^ all energies 
and time scales (particularly the RD times) are given in 
units of Tifio and Hq ^ respectively. The parameters af- 
fecting the numerical results are, v, A, T for the thermal 
noise, a for the system-noise bare coupling and the dipole 
matrix elements {z)nm for the pure MLS. The noise spec- 
trum is assumed to be continuous of which the real part 
is responsible for RD effects and is given by the spec- 
tral density in Q . In the Markovian limit the imaginary 
part of T{t — t') is vanishingly small and the resulting 
Lamb-type energy renormalization effects are negligible. 
In the numerical calculations however we include the full 
complex noise correlation function asii 



10" 



10"' 



■ 10' 



10"' 



10 





• p„ 


(v^-l.cc^.OI) 






\pj 


(v=-l. 0=0.01) 






— Pn 


(v^O, crf.Ol) 






IP,ol 


(v-0, ot-0.01) 






• P„ 


(v^-l.o^.l) 






Ip.oI 


(v=-l. 0=0.1) 






p„ 


(v=-l. 0=0.5) 






lp,.l 


(v=- 1.0=0.5) 





20 



40 60 
time X SI,, 



80 



100 



FIG. 5: Time dependence of the RDM in units of Qq^ for 
various representative a, v parameter sets at T = and A = 
0.1. For the model system the potential is symmetric and the 
bare TLS is in SD configuration. 



X |coth(Lj/2T) cosLj{t -t')~i sin uj{t - i')} 

(11) 

Inserting Hll|l in Q we obtain the system-noise kernel 
for our model. An upper frequency cutoff of LOmax = 5A 
is used in the numerical integral in Hll|) . 

B. Overview of the 2LS results 

We now examine the time behaviour of the RDM in 
© for a 2LS. The solution is shown in Fig[Slon the log- 
arithmic scale for a few representative parameters and 
for the SD configuration. The degeneracy parameter is 
r]2 ~ 10^. We also fixed a ~ 0.01 in the rest of the work 
unless otherwise stated. 

In Fig. the first observation is that, exponential 
RD is effective immediately in the short time regime 
t < 2051q ^. We also confirmed numerically that the 
asympthotic time behaviour as well as the RD rates 
are independent of the initially prepared state. As the 
asymptotic time behaviour is concerned, for symmet- 
ric configurations (pure coupling), the density ma- 
trix converges to the maximum entropy (informationless) 
limit 1/2, where / is the unit matrix, irrespectively from 
the spectral properties of the noise or the system-noise 
coupling. The results also indicate that the relaxation 
time scale (read from the filled symbols) and the de- 
phasing time scale Tdep (read from the hollow symbols) 
are compatible. This result is in agreement particularly 
with the recent exact 2LS calculations using the path 
integral influence functional formalism'*. 

We also confirm that all regions in the noise spec- 
trum have strong influence on RD. For this observation. 
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time X 

FIG. 6: Time dependence of the RDM when the system po- 
tential is biased by ^ptias ~ 0.2. The other parameters are 
a = 0.01, A = 0.1 and u — —1. Comparing this figure with 
the previous one, a crossover in the time dependence of the 
RDM is observed in the u — T plane. The v = Q data is also 
added in the inlet to facilitate this comparison. 



one has to compare Fig's. I5I9I correponding to different 
spectral properties and temperatures. For instance, for 
a = 0.01, A = 0.1 and T = [see Fig|S], we recover the 
under damped and weak dephasing limit of'^ for all i/. In 
addition, no oscillations are observed in the SD configu- 
ration. 

Larger RD rates are observed as i' is made to be more 
negative towards i/ = —1. We identify this behaviour as 
the manifestation of the region-II in the noise spectrum 
(see the end of section II. A). More data are also shown 
in the same figure for indicating the effect of various a 
values. In Figl^the nondegenerate case with an energy 
difference between the levels AE = 1.27iOo is shown for 
a — 0.01,1^ = — 1, A = 0.1 and for various temperatures. 
The V = and v = I curves again yield weaker RD rates 
within the indicated temperature ranges. Another noti- 
cable feature in Figgis the weakly oscillating behaviour. 
The weak oscillations are more prominent at high tem- 
peratures and short times. In this case for lu ^ T the 
system-noise coupling is large due to the large number 
of thermally activated environmental modes. This be- 
haviour, which is characterized by weakly damped Rabi- 
like oscillations, was predicted in the analytic calcula- 
tions of Leggett et al. In order to examine the influence 
of the spectral width A a similar zero temperature plot 
as in Fig[Sl is made in Fig|7| for a wider spectral width 
A = 1. Comparing this figure with Fig[Sl a crossover 
in the time dependence of the RDM can be observed in 
the A — 1/ plane. (The crossover can also be activated 
thermally as to be seen in the following Fig.'s|Hland|^) 

There are a number of ways to increase the effective 
system- noise coupling. The following Fig|Sl gives a sam- 



FIG. 7: Time dependence of the RDM in the SD configuta- 
tion. The other parameters are fixed at T = 0, a = 0.01 and 
A = 1. 



pie from the a — T behaviour for the asymmetric poten- 
tial configuration. Increasing the temperature increases 
the coupling by filling the available photon modes for 
u! < T. The sample data is shown for T = 0,1,5 at 
a = 0.01, A = 5 and = — 1 (indicated by the filled and 
opaque symbols connected with solid lines). The second 
way to increase the system-noise coupling is to directly 
increase a (indicated by the dotted dashed lines). In 
the small coupling regime, for which the sample data is 
shown for a = 5 x lO"'', 10"^^ 5 x 10"^ at T = = -1 
and for an increased A (A = 5), the RDM experiences 
stronger damping. The weakly damped oscillations sur- 
vive at short times at finite and small temperatures. The 
larger the temperature the larger the amplitude of the os- 
cillations and the faster they diminish. 

For completeness we also add in FigElthe behaviour 
in the A — v plane at zero temperature. A comparison 
between the Fig's |S1 and [S] reveals a temperature mod- 
ulated crossover (confirm, for instance, a similar decay 
of the sets for A = 5, = -1 at T = 5 in Fig |H1 with 
A = 5, = 1 at T = in FigEJ. 

The major difference of the model interaction Hamilto- 
nian in from the standard (cr^-type) spin-boson model 
is in the manipulation of the potential. In contrast to 
the standard spin-boson model, in our case only non- 
diagonal, ax, type coupling is present under the sym- 
metric potential [see FigEj. As a result, dramatic dif- 
ferences in the time dependence of the reduced system 
are observed between the two models. For instance, the 
diagonal coupling is standardly considered for the study 
of pure dephasing. In this type of coupling the relax- 
ation is manifestly forbidden and the initial states do 
not change their populations. The diagonal coupling 
also yields strongly temperature dependendent dephas- 
ing rates with the rates vanishing at T = 0. On the 
other hand, when the system-noise coupling is not diag- 
onal, the induced transitions between the system states 
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FIG. 8: The a - T behaviour of the RDM is shown for an 
asymmetric potential with iptias = 0.2. The noise parameters 
are a = 10~^, A = 5 and u = —1. 
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pendence of RD time scales on the spectral width A. A 
naive expectation is that for AE ^ A and at very small 
temperatures the resonant transitions are unfavoured and 
there are no environmental modes available therefore the 
relaxation should be inhibited. The point that is of- 
ten missed in this popular argument is the different role 
played by the short time nonresonant transitions. The 
resonant transitions are favoured when the system in- 
teracts with the environment at sufficiently large times. 
The system however relaxes differently at short times by 
prefering to stay off-resonant in its interaction with the 
noise field thereby sampling all regions of the noise spec- 
trum. This causes the strong dependence on A we observe 
at short times. 

In summary, it is confirmed that the rich transient ef- 
fects are observed usually in the short time behaviour in 
which all energy scales in the noise spectrum take part 
rendering the relaxation process sensitive to the relative 
magnitudes of those scales. We confirmed the several 
crossover regions that have been predicted in the path 
integral influence functional calculations. 

The decoherence and dcphasing dynamics is governed 
by all frequency regions in the noise spectrum. In partic- 
ular, the short time behaviour is affected strongly by all 
frequency regions due to the nonresonant processes. For 
A ^ 1, Fig|Sl indicates that as v increases in the interval 
— 1 < < f the relaxation rates gradually increase still 
remaining in the weak relaxation regime. A crossover in 
the v — A plane is observed [compare with Fig[7] as A is 
increased. With this being the case for symmetric poten- 
tials, for asymmetric ones the additional feature of weak 
oscillations are present in the short time dynamics. 
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FIG. 9: The density matrix parameterized by different A and 
1/ values at T = and a — 0.01 for the asymmetric potential 
with ifbias ~ 0.2. The crossover from the weakly oscillat- 
ing damped behaviour at short times to strong relaxation at 
longer times can be observed. 



can probe the entire noise spectrum creating decoherence 
even at zero temperature. These induced transitions are 
nonresonant and they have observable effects particularly 
in the short time dynamics of the RDM^^. The RD times 
observed as the result of such system-noise coupling are 
expected to be nonzero even at zero temperature. This 
characteristic behaviour of the non-diagonal coupling is 
confirmed in our calculations both for the 2LS in FigsIS]- 
|51 and for the MLS in the following sections. Recently, 
there are other claims using realistic models on deco- 
herence effects in mesoscopic systemsiS as well as a few 
experimental confirmations on the saturation of the RD 
rates at low temperatures^". 

A curious observation in Fig's I5I9I is the strong de- 



C. MLS with 3<N 

The effect of multileveledness on decoherence has not 
yet received the attention that it deserves in the litera- 
ture. This is, in part, due to the lack of practical analytic 
tools in solving the master equation. The complexity 
of the formal methods such as the noninteracting blip 
approximation increases at each time step as N'^ which 
renders the analytic sum over all virtual configurations in 
the path integral approach intractable. Usually, the com- 
mon argument is that, for sufficiently low temperatures, a 
multilevel system, of which the first two levels (the qubit) 
are sufficiently well separated from the rest, behaves as 
a two state system. We have already observed that there 
are two major pitfalls in this assumption. Firstly, it ex- 
cludes the very realistic case in which decohering effects 
are induced through interactions with a strongly fluctuat- 
ing quantum field. In such a case, the fluctuations in the 
distribution of environmental modes in the noise spec- 
trum is the major source of decoherence. Secondly, and 
more generally, the short time behaviour of a MLS -which 
is the most prominent regime in the quantum computa- 
tional perspective- is affected by a large frequency region 
in the noise spectrum. These comprise the basic motiva- 



8 
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FIG. 10: Comparison of the eflect of the spectral width A on 
the time dependence of the elements poo , pii and |poi | between 
the 2LS and 3LS. The fixed parameters are T = 0, a = 0.01 
and v = —1. The open symbols refer to the case N = 2 
and the solids ones refer to A'^ = 3. Also the solid lines are 
the diagonal elements poo (with poo(O) = 0.1) and pn (with 
Pii(O) — 0.9), the dotted-dashed lines are the non-diagonal 
ones IpioI (with |pio(0)| = 0.3). More specifically, circle is 
A = 1, square is A = 10. The inlet is the case v = 1 and only 
A = 10 is displayed for simplicity. 

tion of why we should look into the effect of higher levels 
on decoherence in a MLS. We will continue by examining 
the cases = 3 and 3 < iV separately. In particular, we 
remark on certain interesting decoherence properties of 
the doubly degenerate systems with 4 < A^. 

A^ = 3 case 

From the quantum computational point of view, the 
three level systems are as important as the two level ones 
(see for instance^^). Recently, the Linblad approach was 
used for the RDM of a multilevel system in connection 
with the quantum Zeno effectS-^ . The Linblad equation is 
derived directly from the Bloch-Redfield-Fano equation 
for the RDM and is based on the validity of the Marko- 
vian condition A"-'^ <C tr. We compare in Fig^| the 
solutions of the RDM for A'' = 2 and A^ = 3 at a = 0.01, 
T = and for A = 1, 10. We look at the symmetric po- 
tential in the SD configuration. For 3LS, the observed 
energies form a A-shaped configuration and, in units of 
^0, are roughly Eq ~ Ei k. 0.1 and E2 « 2.3. The 
degeneracy parameter 77 = {E2 — Ei)/{Ei — Eq) ~ 10^ 
and the third level E2 is above the double well barrier as 
shown in Fig. Here we mainly discuss the subohmic 
case v = —I. 

When A <C 1 the resonant coupling of the first two lev- 
els to the third level is very weak. In principle, at short 
observational times, the nonresonant excitations are in- 
duced with frequencies much higher than the resonant 
frequency E3 — E2 — 2 (in units of Qq). When A ^ 1 



however, these transitions are suppressed by the Gaus- 
sian cutoff. As a result, the 3LS is basically confined to 
its highly degenerate qubit subspace. This confinement 
can be observed all the way up to much higher spectral 
widths such as A = 1 as long as the Gaussian suppres- 
sion is manifested. This behaviour is shown in Fig llOl for 
A = 1 and A = 10. For considerably long duration (i.e. 
~ 100 X ^Iq ^) and for A = 1 the qubit subspace in the 
three level system has a negligible leakage into the third 
level. For much larger A such as A = 10, the third level 
is allowed to participate in the transitions. The RD rates 
are therefore found to be significantly larger than that of 
the N = 2 case before. 

An enhancement in the high frequency and suppression 
in the low frequency coupling as compared to = — 1 , is 
generated if we now consider j/ = 1 (shown in the inlet 
of Fig llO|) . Under these conditions, the short time be- 
haviour is dramatically influenced by a strong leakage to 
the third level (note the short time span on the horizontal 
axis between zero and 5 x f^g ^). 

The asymptotically long time dynamics is indepen- 
dent of the system-noise parameters. In sufficiently long 
time the system looses all the information that is put in 
the initial state: pkkit 00) = 1/3, (fc = 0,1,2) and 
|pfe,(t^oo)| =0,(fc^j). 

A^ = 4 case 

We compare the 4LS with a 2LS in Fig^Jfor v = — 1. 
Here, we have three sets of curves indicated by (a), (b) 
and (c). In Fig^Ja the 4LS is compared to 2LS when 
both systems are highly degenerate. In Fig^Jb the four 
level system is doubly degenerate and at zero temper- 
ature. The third set of curves are plotted in Fiallllc 
corresponding to the DD configuration at finite temper- 
atures. 

For the singly degenerate (SD) case, the qualitative 
features between the 2LS and the 4LS are similar to the 
previously discussed case between 2LS and the 3LS. Here, 
we observe for the diagonal elements, a much higher re- 
laxation rate (as well as leakage) out of the qubit sub- 
space during the observed time although the dephasing 
rates are indistinguishable for the 2LS and the 4LS. The 
RDM asymptotically approaches to the informationless 
limit p — 1/4. 

For A^ = 4 we have the chance to look at the DD 
configuration as depicted in Fig. |j2Jl. For the doubly 
degenerate configuration, the degeneracies are as high 
as ri2 = {E2 - Ei)/Ei - Eq) ~ 3 x 10^ and 771 = 
{E3 — E2)/E2 — El) ~ 10^. We surprisingly observe a 
tremendous suppression (by almost two orders of magni- 
tude) at sufficiently small temperatures in the RD rates 
(Fig^Jb). The rates and the DD-suppression strongly 
depend on the temperature. A comparison of the T = 
0,A = 1 in Figinib and T = 1,A = 1 in FigElc can 
reveal this strong dependence. 

4 < N case 
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FIG. 11: (a) Comparison, at T = and v = —1, of the 
effect of the spectral width A on the time dependence of the 
elements poo,/9ii and |poi| between the 2LS and 4LS. The 
symbols are the same as in Fig llUI The figure describes the 
singly degenerate case r\ ~ 10®; (b) the rates at T = and 
V = —\ for the doubly degenerate (denoted by DD in the 
figure title) configuration r\2 — (£2 — E\)/E± — Eo) ~ 3 x 10® 
and ryi = (iJs — E2) / E2 ~ Ei) 10® [horizontal and vertical 
axes have the same span as (a)]; (c) the thermal case at the 
indicated A, T values at v = —1. 



In order to extract some quantitative numbers for the RD 
times we made use of the numerical observation that for 
weak system-environment coupling the time dependence 
is approximately exponential at short times. We then 
follow'^ and write for the time dependence of a general 
matrix element at short times 

\pij{t)\ ~ |py(oo)| + [|py(0)| - |py(oo)|]ea;p(-i/Ty) . 

(12) 

The RD times are extracted from the time dependence 
of pii(t) and |pio(t)| respectively as 



1 



dh\\pij\ 



1 - \pij{(x,)/ pij{0)\ dt 



(13) 



where i = j = 1 is used in the calculation of the re- 
laxation rate and i = 0, j = 1 is for the qubit dephas- 
ing rate corresponding to the first excited level. For the 
RDM at asympthotic times we have /5ii((X)) = 1/N and 
|pio(oo)| = 0. The equation ifT^ breaks down when 
|Pij(0)| = |Pij(oo)l which we stay away by appropri- 
ately choosing pij{0). In FigElthe data are represented 
at zero temperature and 1^ = 0. To be used in (|13|l . 
the initial conditions are set at jOoo(O) — 0.2,pii(0) = 
0.8,pio(0) = 0.4i with all density matrix elements out- 
side the qubit subspace zero. Three different curves stand 
for (bottom to top) A = 0.1, 1, 10 with the open symbols 
corresponding to dephasing and the solid ones to the re- 
laxation rates. Each set of data is shown for SD as well as 
DD configurations separately. Also note that the vertical 
axis is logarithmic. 



FIG. 12: Relaxation and dephasing rates against the number 
of levels for different spectral widths at T = and = 0. 
Note the logarithmic vertical axis. Small symbols refer to 
the singly degenerate MLS and the larger symbols refer to 
doubly degenerate one. The open and solid symbols refer to 
dephasing and relaxation times respectively. 



Let us concentrate first on the SD configurations in 
FigEl In a large A range N = 4 and ~ 10 appear 
to be two crucial points. For 4 < iV relaxation is ap- 
proximately twice faster than dephasing and both rates 
rapidly saturate near A'^ ~ 10 and they are independent 
of A^ for 10 < A^. The onset of saturation is naturally 
model dependent. In our case this onset coincides with 
the range of strong dipole transition matrix elements of 
the model in (QJ (see Fig|3J). Turning to the DD con- 
figurations, we observe that for the same environmental 
parameters and for all A^, decoherence rates for the DD 
case are strongly suppressed by nearly two orders of mag- 
nitude as compared to the SD configuration. 

With this section we conclude the RD calculations for 
the interaction between the system and the thermally 
equilibrated noise. We now focus our attention on the 
calculations of the transition rates by a different ap- 
proach, the Fermi Golden rule. 



IV. FERMI GOLDEN RULE 

The Fermi Golden Rule (FGR) provides a simple and 
qualitative tool to reproduce many of the features of the 
relaxation times that we observe in Fig ll2l Quantitative 
agreement should not be expected between the Figll2l 
and the FGR results. This is mainly due to the fact 
that the data produced in Fig^l reflects the effects of 
the short time dynamics whereas, the FGR gives more 
accurate results for the long time resonant interactions. 
In comparing the time scales found by directly solving 
the RDM and by FGR, we keep the absolute time scales 
arbitrary and only compare the qualitative behaviour. 
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FIG. 13: Relaxation rates against the number of levels for the 
SD (solid symbols) and DD (open symbols) cases. 



We assume that the MLS is prepared at t = in the 
first excited state |{C}, !)• The probability that the sys- 
tem stays in the same state after interacting with the 
environment for a duration t is- 



FIG. 14: Comparison of the absolute relaxation rates (in com- 
mon arbitrary units) with A'^ between the SD (solid symbols) 
and DD (open symbols) configurations. Here, u — and 
A = 10. Same temperatures are implied for each SD-DD 
pair. 



dt'n^nt{t')]\ij{0)} (14) 



where for our case \ip{0)) = |{C},1). Including second 
order perturbation in the dipole couplings with an en- 
vironment in thermal equilibrium, 1)14(1 can be written 



1 -P4>it) 



dt' 

^0 



dt"e'(^"o-^^=)(t'-*")^(t'_t") 

(15) 



where !F{t' — t") is the environmental correlation function 
given by (|10|l . The relaxation rate t^^j^ corresponding 
to the first excited state is found by the same method 
that is described in {^) [The Eq. (3.37) therein]. 

The relaxation times found by the FGR are summa- 
rized in Fig' dl3l and 1141 In Fig^|the relaxation times 
are plotted against the same parameters as in Fig ll2l with 
the same symbols. The FGR data reproduce many of the 
features in FigEl The first observation is the same offset 
at = 4 and the saturation of the rates slightly above 
this offset with a rapid increase for 4 < iV. The second 
observation is that by increasing the spectral width, the 
relaxation rates can be increased by as much as an order 
of magnitude. 

When the temperature is varied, we observe the same 
trend as in the previous figure as depicted in Fig. I|14|l . 



CONCLUSIONS 



We examined the RD properties of a multilevel sys- 
tem using non-Markovian master equation formalism. It 
is shown that the short time behaviour of the density 
matrix is influenced by nonresonant transitions in the 
MLS receiving contributions from all frequency regions 
in the noise spectrum. For the model interaction used, 
the dipole transitions are nonzero within a finite range of 
levels. The RD times calculated within this model show 
a saturation within the same range largely independent 
from the system- noise parameters. 

It is generally found that the decoherence effects in 
MLS are more pronounced than those in the 2LS. We 
observe that a distinct counter example is posed by the 
doubly degenerate MLS with 4 < N. The RD rates are 
found to be highly suppressed in comparison with the 
singly degenerate or non-degenerate systems for the same 
system and environment parameters. These result were 
also confirmed using the transition rates found from the 
Fermi-Golden rule. At the first glance, this curious sup- 
pression of decoherence reminds us the decoherence free 
subspaces. Nevertheless, the arbitrariness of the param- 
eters of the DD model, and in particular of the dipole 
matrix elements rules out the possibility whether any set 
of invariant states can form a decoherence free subspace 
under the coupling with the environment. In order to un- 
derstand the true nature of this strong suppression, more 
formal and analytic methods must be developed for the 
DD systemsS^. 
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